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\ We describe a simple method of determining whether the singularity that forms in the 

spherically symmetric collapse of inhomogeneous dust is naked or covered. This derivation 

considerably simplifies the analysis given in the earlier literature, while giving the same 

0\ , results as have been obtained before. 
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Various authors have shown that the spherical gravitational collapse of inhomogeneous 
dust results in the formation of a curvature singularity which is naked for certain initial 
conditions, and covered for other initial conditions. Probably the first (numerical) results 
on this problem are due to Eardley and Smarr |l[ . This was followed by the analytical 
work of Christodoulou and Newman ||] who considered smooth initial data. Their results 
were generalised, among others, by Ori and Piran and by Dwivedi, Jhingan, Joshi and 
Singh H, ||. As is only to be expected, in one way or the other, these works all deal with 
propagation of null geodesies in the spacetime of collapsing dust. 

While successive works have succeeded in simplifying the earlier analysis, perhaps it can 
be said that the discussions continue to remain somewhat involved. In the present paper, 
we describe a short but straightforward method of showing whether the naked singularity 
in the dust model is covered or naked. We reproduce results obtained previously by other 
methods. We restrict attention to the marginally bound dust collapse - similar principles 
may be used to derive results for the non-marginally bound case. 

In comoving coordinates (t, r, 8, (f>) the spacetime metric for spherical dust collapse is 
given by 

ds 2 = dt 2 - R ,2 dr 2 - R 2 dQ 2 (1) 

where R(t, r) is the area radius at time t of the shell having the comoving coordinate r. A 
prime denotes partial derivative w.r.t. r. The energy-momentum tensor for dust has only 
one non-zero component Tq = e(t,r), which is the energy density. The Einstein equations 
for the collapsing cloud are 

c 4 e[t,r) Rm>> K r ■ { ) 

A dot denotes partial derivative w.r.t. time t. The function F(r) results from the integration 
of the second order equations. Henceforth we shall set 87rG/c A = 1. 

The second of these equations can be easily solved to get 

R V^ t ,r)=r 3 / 2 --VFt (3) 
2 
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where we have used the freedom in the scaling of the comoving coordinate r to set -R(0, r) = r 
at the starting epoch of collapse, t = 0. It follows from the first equation in @ that the 
function F(r) gets fixed once the initial density distribution e(0, r) = p(r) is given, i.e. 

F(r) = J p(r)r 2 dr. (4) 

Hence F(r) has the interpretation of being twice the mass to the interior of the shell labeled 
r. If the initial density p(r) has a series expansion 

P(r) = po + pir + ^p 2 r 2 + ^p 3 r 3 + ... (5) 

near the center r = 0, the resulting series expansion for the mass function F(r) is 

F(r) = F r 3 + F ir 4 + F 2 r 5 + F 3 r 6 + ... (6) 

where F q = p q /q\(q + 3), and q — 0, 1, 2, 3.. We note that we could set p\ — without in any 
way affecting the conclusions of this paper. Further, the first non-vanishing derivative in the 
series expansion in (||) should be negative, as we will consider only density functions which 
decrease as one moves out from the center. 

According to (|3|) the area radius of the shell r shrinks to zero at the time t c (r) given by 

2r 3/2 

Ur) = -j—. (7) 

At t — t c (r) the Kretschmann scalar 

F' 2 F F' F 2 

diverges at the shell labeled r and hence this represents the formation of a curvature sin- 
gularity at r. In particular, the central singularity, i.e. the one at r = 0, forms at the 
time 

2 2 

to = T7 ^ = ~yr=- (9) 

At t — to the Kretschmann scalar diverges at r = 0. Near r = 0, we can expand F(r) and 
approximately write for the singularity curve 

t c (r) = t - -^r». (10) 
3-^0 



Here, F n is the first non- vanishing term beyond Fq in the expansion (|^). We note that 
t c (r) > to, since F n is negative. 



We wish to investigate if the singularity at t = to, r = is naked, i.e. are there one 
or more outgoing null geodesies which terminate in the past at the central singularity. We 
restrict attention to radial null geodesies. Let us start by assuming that one or more such 
geodesies exist, and then checking if this assumption is correct. Let us take the geodesic to 
have the form 

t = t + ar a (11) 

to leading order, in the t — r plane, where a > 0, a > 0. In order for this geodesic to lie in 
the spacetime, we conclude by comparing with fllPP that a > n, and in addition, if a = n, 
then a < -F n /3F 3/2 . 

As is evident from the form ([!]) of the metric, an outgoing null geodesic must satisfy the 
equation 

(12) 



dr 



In order to calculate R' near r = we first write the solution (|3|) with only the leading term 
F n retained in F(r) in This gives 



/? r[l--\lF 



F 
2F n 



2/3 



(13) 



Differentiating this w.r.t. r gives 



F 

2F . 



/i" V3 (i-lMt-^ F ^n. 



A\/F< 



(14) 



o 



Along the assumed geodesic, t is given by (|TT|). Substituting this in R' and equating the 
resulting R' to dt/dr = aar a ~ l gives 



aar 



a-l 



1 - IVK [t + ar a ] - ^5^r« [t + ar< 



1-IVFo 



1 + 



[t + ar c 



1/3 



(15) 



This is the key equation. If it admits a self-consistent solution then the singularity will be 
naked (i.e. at least one outgoing null geodesic will terminate at the singularity), otherwise 



not. We simplify this equation by putting in the requirement mentioned earlier, that a > n. 
Consider first a > n. In this case we get, to leading order 

aar^=(l + ^)(-^-) 2/ \^ (16) 



3 ) V 2F , 

which implies that a = l+2n/3, and a = (—F n /2F ) 2 ^ 3 . By substituting integral values for n 
we find that only for n = 1 and n = 2 the condition a > n is satisfied. Hence the singularity 
is naked for n = 1 and n = 2, i.e. for the models p\ < and for p\ = 0, p 2 < 0. There is at 
least one outgoing geodesic given by ( |TT| ) which terminates in the central singularity in the 
past. If n > 3 then the condition a > n cannot be satisfied and the singularity is not naked. 
This is the case p\ = p 2 = P3 = 0. 



nar n - 1 = 1 _ - ^ r 2 "/ 3 (17) 



Consider next that a = n. In this case we get from (|15|) that 

3„ rW (2»+3)F n 
( F n 3a rnr 

which implies that n = 3 and gives an implicit expression for a in terms of F 3 and -F . This 
expression for a can be simplified to get the following quartic for a: 

12/^a 4 - a 3 (-4F 3 /F + F 3/2 ) - 3F 3 a 2 - 3F 3 2 /F 3/2 a - (F 3 /F ) 3 = 0. (18) 

By defining b = a/F and £ = F 3 /F^ 2 this quartic can be written as 

46 3 (36 + - (6 + 3 = 0. (19) 

The singularity will be naked if this equation admits one or more positive roots for b which 
satisfy the constraint b < — £/3. This last inequality is the same as the condition a < 
—F n /3F^ 2 given below equation fllTf ). We note that £ is negative. This quartic can be made 
amenable to further analysis by substituting Y = — 26/£, and then r\ = — 1/6£, so as to get 

Y'\Y - 2/3) - 7](Y - 2) 3 = 0. (20) 

As discussed in || this quartic has two positive real roots provided 77 > T]i or 77 < 772 where 

26 26 
7& = - + 5 VS, 772 = J - 5^3. (21) 

5 



We also require that Y < 2/3. By examining the quartic (B^) one can see that if r\ > r\\ then 
Y > 2; hence this range of rj is ruled out. Thus the singularity is naked provided 77 < r] 2 , or 
equivalently £ < -25.9904. 

This completes the analysis to decide whether or not the central singularity is naked, and 
we get the same results as have been given earlier in the literature, albeit in a much simpler 
manner. Now we examine whether or not there is an entire family of radial null geodesies 
which terminate at the naked singularity. For this purpose we assume a solution for the 
geodesies correct to one order beyond the solution (|TT|), i.e. we take 

(22) 



t = t Q + ar a + dr a+l3 . 



where d and (3 are constants to be determined, and a and a take the values calculated above. 
As before, we substitute this form of t{r) in the expression (|i~4T) for R' to get 



R 



2VF0 



t + ar a + dr a+ P 



(2n+3)F n , 
4^ ' 



to + ar a + dr a+p 



(1 



En 
2F ' 



1 + TTFrr' 



[t + ar a + dr 



o+/3 



1/3 



(23) 



Next, we equate this R' to the dt/dr calculated from (p2|). For the cases n = 1,2 we get, 
after retaining terms up to second order 



aar 



a-l 



+ (a + (3)dr^- 1 =(l + ^ 



2Fn 



2/3 



„2n/3 1 f) r a - n /3^ 



(24) 



As before, at the leading order a and a get fixed. At the next order, we get j3 = l + n/3 and 
d = D/(2 + n/3) where 



D = -JF 
2 V 



2F 



1/3 



1 / 2n 
■1 + - 1 + — 
3 V 3 



2F 



-2/3' 



(25) 



It thus follows, according to (0), that when n — 1, 2 there is to this order only one outgoing 
geodesic, having the values of d and f3 given above. 



Consider next n = 3. By repeating the above calculation of R' we get 

. r 2+/3 + o(r 5 ). 



R' = h^F e + \ /3 r 2 + 3.2^bd 

2 (f + 36) 1/3 {i + 3bf 3 



(26) 



Here, 0(r 5 ) is a term of order r 5 which is independent of d. Further analysis depends on 
whether or not (3 is less than 3. Assume first that (3 < 3. Then the 0(r 5 ) term can be 
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ignored. By equating R' to dt/dr = 3ar 2 + (3 + j3)dr 2+/3 , we get the earlier quartic fllTf ) for 
b. At the next order, d drops out and one gets an equation for j3, i.e. 

3 + /? = 3&.2 1 / 3 (27) 

(£ + 36) 4/3 

Since d drops out, this means that d is arbitrary, and there will be an entire family 
of outgoing null geodesies terminating at the singularity, provided (3 is non-negative. It is 
essential that (3 be non-negative, otherwise these geodesies will not lie in the spacetime, as 
is evident from a comparison with the singularity curve (|T0|) . As we saw above, the quartic 
( pUD and hence ([19|) has two positive roots in the naked singular range. We now show that 
(3 as defined in (|27| ) is positive at one of the roots, and negative at the other root. Let us 
write the quartic ( f20|) as V(Y) = where 

V(Y) = Y 3 (I - y) - a(2 - F) 3 . (28) 



^3 

It is then easily shown that 



'dV\ _,o/2 
dY 



PY 2 (~ - n) (29) 



Y=Y 

where Y Q is a real, positive root of the quartic. Since the derivative dV/dY must be positive 
at one of the roots and negative at the other, and since Y < 2/3 it follows that j3 is positive 
at one of the roots and negative at the other. Hence one of the roots admits only one 
outgoing geodesic (for which d=0) while the other root admits an entire family of outgoing 
geodesies. It is easily verified that the family emerges from the larger of the two roots, which 
lies closer to the singularity curve (p~0|) . 

However, it also turns out, as is verified numerically, that the positive value of f3 does 
not remain below 3 for all £. It can be shown that (3 can be written as 

Y 2 -8Y + A 

fj ~~(2/3-Y) (2-YY m 

For values of £ smaller than a certain critical value, j3 becomes larger than 3, so that then 
the 0(r 5 ) term in fl2"E|) dominates over the term proportional to r 2+/3 . For such cases, we get 



from dt/dr = R' that (5 — 3, and d also gets fixed at a particular value. In order to see the 
family of outgoing rays we will have to look at higher order terms in the various expansions. 

Consider next the case of ingoing rays, given by dt/dr = —R', for which we take 

t = t - er 3 - gr 3+1 (31) 

(We consider only the case n — 3). The expression for R 1 is 



R' = V 3 F - e i/3 r 2 + ^2 4 / 3 l -^r^ (32) 



where h = eF . Equating this to dt/dr = — 3er 2 — (3 + 7)gr 2+7 gives, at order r 2 , the 
following quartic for h 

h 3 (12h - 40 + (h - 3 = (33) 

which admits a positive root /i for all £, as expected, so that there is always an ingoing ray. 
However at the next order, we get the relation 

-(3 + 7) = ^2 4 / 3 (34) 

4 (£-3h) 4/3 

which cannot be satisfied, unless g — 0, since the l.h.s. is negative and the r.h.s. positive. 
This shows that there is only one ingoing ray to all orders in the expansion. 
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